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The interaction of compact objects with an infinitely extended mirror plane due to quantum 
fluctuations of a scalar or electromagnetic field that scatters off the objects is studied. The mirror 
plane is assumed to obey either Dirichlet or Neumann boundary conditions or to be perfectly 
reflecting. Using the method of images, we generalize a recently developed approach for compact 
objects in unbounded space [H [2] to show that the Casimir interaction between the objects and 
the mirror plane can be accurately obtained over a wide range of separations in terms of charge 
and current fluctuations of the objects and their images. Our general result for the interaction 
depends only on the scattering matrices of the compact objects. It applies to scalar fields with 
arbitrary boundary conditions and to the electromagnetic field coupled to dielectric objects. For 
the experimentally important electromagnetic Casimir interaction between a perfectly conducting 
sphere and a plane mirror we present the first results that apply at all separations. We obtain 
both an asymptotic large distance expansion and the two lowest order correction terms to the 
proximity force approximation. The asymptotic Casimir-Polder potential for an atom and a mirror 
is generalized to describe the interaction between a dielectric sphere and a mirror, involving higher 
order multipole polarizabilities that are important at sub-asymptotic distances. 



I. INTRODUCTION 

The interaction between neutral objects is dominated 
by fluctuation forces due to the coordinated behavior of 
fluctuating charges or collective modes inside the objects. 
At zero temperature or at sufficiently small distances, the 
interactions result from quantum fluctuations. They are 
important on atomic scales as well as between macro- 
scopic bodies. A prominent example for the latter is 
the Casimir force between parallel metallic plates due to 
current or electromagnetic field fluctuations [3J. During 
the last decade experimental verifications of this effect 
have been performed with increasing precision. These 
high precision measurements were enabled by the use of 
curved surfaces instead of parallel planar mirrors in order 
to avoid the problem of parallelism. The most commonly 
employed geometry is a sphere-plate setup that was used 
in the first high-precision tests of the Casimir effect [4] [5j . 
This geometry has been successfully used ever since in 
most of the recent experimental studies of Casimir forces 
between metallic surfaces H El 13 EH HI1 EH HH]- In 
order to keep the deviations from two parallel plates suffi- 
ciently small, spheres with a radius much larger than the 
surface distance have been used. The effect of curvature 
has been accounted for by the "proximity force approxi- 
mation" (PFA) [14 . This scheme is assumed to describe 
the interaction for sufficiently small ratios of radius of 
curvature to distance. However, this an uncontrolled as- 
sumption since PFA becomes exact only for infinitesimal 
separations, and corrections to PFA are generally un- 
known. 

At the other extreme, the interaction between a pla- 
nar surface and objects that are either very small or at 
asymptotically large distance is governed by the Casimir- 
Polder potential that was derived for the case of an atom 



and a perfectly conducting plane [IS] . This limit has been 
probed experimentally with high precision for a Bose- 
Einstein condensate that was trapped close to a planar 
surface [IB] . Recently, there have been first attempts to 
go beyond the two extreme limits of asymptotically large 
and small separations by measuring the Casimir force be- 
tween a sphere and a plane over a larger range of ratios 
of sphere radius to distance [17] . 

So far, no theoretical prediction is available that can 
describe the electromagnetic Casimir interaction between 
a compact object and a planar surface at all distances, 
including the important sphere-plate geometry. Until re- 
cently, progress in understanding the geometry depen- 
dence of fluctuation forces was hampered by the lack of 
practical methods that are applicable at all separations. 
Conceptually, the effect of geometry and shape is dif- 
ficult to study due to the non-additivity of fluctuation 
forces. Explicit consequences of this non-additivity and 
also non-monotonic changes in the force have been re- 
cently predicted for a pair of cylinders next to planar 
walls [TB]- This behavior has been interpreted in terms 
of collective charge fluctuations inside the bodies 1 '.) . 

For some decades, there has been considerable inter- 
est in the theory of Casimir forces between objects with 
curved surfaces. Two types of approaches have been pur- 
sued. Attempts to compute the force explicitly in partic- 
ular geometries and efforts to develop a general frame- 
work which yields the interaction in terms of characteris- 
tics of the objects like polarizability or curvature. Within 
the second type of approach, Balian and Duplantier stud- 
ied the electromagnetic Casimir interaction between com- 
pact perfect metals in terms of a multiple reflection ex- 
pansion and derived also explicit results to leading order 
at asymptotically large separations [20] . For parallel and 
partially transmitting plates a connection to scattering 
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theory has been established which yields the Casimir in- 
teraction of the plates as a determinant of a diagonal 
matrix of reflection amplitudes [21] . For non-planar, de- 
formed plates, a general representation of the Casimir 
energy as a functional determinant of a matrix that de- 
scribes reflections at the surfaces and free propagation 
between them has been developed in Ref. [52]. Later on, 
an equivalent representation has been applied to pertur- 
bative computations in the case of rough and corrugated 
plates with finite conductivity [23, 24 . 

Functional determinant formulas have been used also 
for open geometries that do not fall into the class of par- 
allel plates with deformations. For the electromagnetic 
Casimir interaction between a planar plate and infinitely 
long cylinders, a partial wave expansion of the functional 
determinant has been developed [25]- The same results 
have been reobtained and used to compute corrections 
to the PFA for the cylinder-plate geometry in Ref. [25] . 
Kenneth and Klich identified the inverted Green's func- 
tion in the functional determinant as a T-matrix and de- 
rived a formal result for the Casimir interaction in terms 
of this matrix in the case of scalar fields in a medium with 
a space and frequency dependent speed of light [27] . 

Recently we described a new method based on a multi- 
pole expansion of fluctuating charges that makes possible 
accurate and efficient calculations of Casimir forces and 
torques between any number of compact objects [U [2]. 
The method applies to electromagnetic fields and di- 
electrics as well as perfect conductors. It also applies to 
other fields, such as scalar and Dirac, and to any bound- 
ary conditions. In this approach, the Casimir energy is 
given in terms of the fluctuating field's scattering am- 
plitudes from the individual objects, which encode the 
effects of the shape and boundary conditions. An equiv- 
alent partial wave expansion has been applied to a scalar 
model for dielectrics |28j . 

Casimir interactions due to scalar field fluctuations 
serve as a simplified model for the full electromagnetic 
interaction. This model is usually easier to analyze and 
provides an important tool in developing conceptually 
new approaches and in estimating geometry dependen- 
cies. For the sphere-plate geometry with Dirichlet bound- 
ary conditions Bulgac et al. obtained the Casimir inter- 
action over a wide range of separations from a modified 
Krein trace formula [29 , 30 . Most notably in the context 
of scalar fields, a versatile numerical world- line algorithm 
based on Monte Carlo methods has been developed and 
applied to a number of interesting geometries, including 
the here studied sphere-plate interaction [31] 132] . 

In this work, we extend our approach developed in 
Refs. [T][2] to describe the interaction of compact objects 
in the presence of a plane mirror. Our general result 
holds for scalar and electromagnetic fields. In the elec- 
tromagnetic case the mirror is assumed to be perfectly 
conducting but the compact objects can be dielectrics or 
perfect conductors. The derivation of the general result 
for the Casimir interaction, see Eqs. (11.48 1 and (11.100 1, 



gral approach and the method of images. In Section III 
we apply our approach to compute the interaction of a 
sphere with a plane mirror over a wide range of separa- 
tions for the scalar Dirichlet and Neumann problem and 
for the electromagnetic field that is most relevant to ap- 
plications. We provide a large distance expansion of the 
interaction, generalizing the Casimir-Polder potential to 
include higher order multipole polarizabilities. For small 
separations, we compute the two leading correction terms 
to the PFA. 



II. INTERACTION BETWEEN FLUCTUATING 
SOURCES AND THEIR IMAGES 

In this Section we first review the functional integral 
formulation of Casimir interactions between compact ob- 
jects for a scalar field and the electromagnetic field. The 
fluctuating field is integrated out in order to obtain an 
effective action for the fluctuating sources on the ob- 
jects. Then we show that the interaction of compact 
objects with an infinite plane mirror can be described by 
the equivalent problem of the interaction of compact ob- 
jects with their mirror images in otherwise empty space, 
i.e., without the plane mirror. This equivalence holds for 
Dirichlet or Neumann boundary conditions (at the mir- 
ror and the objects) for a scalar field and for a perfectly 
conducting mirror plane and arbitrary dielectric objects 
in the case of electromagnetic fluctuations. 



A. Scalar field 

We first consider a real quantum field <j>(x, t) in order 
to introduce the method of images in the path integral 
formulation of Casimir interactions. We assume that the 
space is divided into two half spaces by a mirror plane 
at z = 0. N fixed closed surfaces S Q , a — 1, . . . , N, are 
located in the right half space (z > 0). These surfaces 
are regarded as the boundaries of objects on which cither 
Dirichlet or Neumann boundary conditions are imposed. 
The action of the unconstrained field in Minkowski space 
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dt / dx 
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(d t <}>) 2 - (v# 



(ii.i) 



Here the x-integration runs over the right half space, in- 
dicated by j > . The free energy T of the constraint field 
</>(x, t) at inverse temperature (3 is represented by the 
Euclidean functional integral 



[v^ce^i-SE^yn) 



{11.2) 



where the Euclidean action Se follows from Eq. (II.I I 



after the Wick rotation ct 
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is given in Section [n] by combining a functional inte- 



SE = fcJ dX J dx {( 9 o0) 2 + W) 2 } (IL3) 
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with A = f3hc. The functional integral is over all fields 
that are periodic in the Xo-interval from to A and that 
obey the boundary conditions at the surfaces, indicated 
by the subscript C. 

The surfaces are fixed and the boundary conditions 
are time independent. Hence each Fourier component of 
the field with respect to x° obeys the constraints at the 
surfaces separately. To make use of this property, we 
expand 0(x, a; ) as 



(H.4) 



with Matsubara frequencies n n — 2im/A and 0_„(x) = 
0*(x). The functional integral splits into independent 
functional integrals over the <p n (x) and the logarithm of 



Z can be written as the sum 
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OO 
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In / \D(j> n ]c exp 



A 

~2hc 



dx {k^ 



(II.5) 
|V0„| 2 } 



In the following we are interested in the limit of zero 
temperature. Then A — > oo, and the sum over n can be 
replaced by the integral ^ dn and ^„(x) is replaced 
by <f>(x, k). Combining positive and negative k, we get 



InZ = 



with 



A 
2^ 



dn \nZ{n) 



(II.6) 



Z(k) = y"[P0(x,«)]c[^*(x, K )] c exp -Ay dx { K 2|0|2 + | V( ^|2} 



(II.7) 



The ground state energy is obtained from Z as 
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he 
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hiZ 



dn \tiZ{k) 



(II. 



The Casimir energy is the difference between the ground 
state energy of the surfaces at their actual distance and 
the sum of the ground state energies of the separate sur- 
faces which is obtained by removing the surfaces to infi- 
nite separation, 



he 
"2^ 



dn\n 



Z{k) 



(II.9) 



where Z m (n) is the partition function for the surfaces 
with infinite separation. In the following we suppress the 



label k. Next, the constraints at the objects are imple- 
mented by functional (5-functions [33J [33] . For Dirichlet 
boundary conditions, = 0, on the surfaces S Q , the con- 
straint functional integral can be expressed in terms of 
an unconstrained integral by using 



[V<t>}c[V4>*] c 

N 



(11.10) 



Jv<t>V<t>* f[ Jvg a Vg* a exp i j dx {g* a (x)0(x) + c.c. } 



where the (5-functions at each position of the surfaces 
have been written as an integral over a source field g a (x) 
that is non-zero on the surfaces E Q only. When we use 
this representation of the constraints in Eq. (II. 7 1, the 
now unconstrained functional integral over (j>(x) is Gaus- 
sian and yields 



Z{k) = Z f[ ( Vg a Vg* a exp [ dx [ dx ' {^( X ) G >( X ' x '' k )pp( x ') + cx - } 



Za II / ' D g a T)g* a t 



-(hc/A)S[g] 



(11.11) 



where Zq is the partition function of the field in the maim boundary conditions at the surfaces S Q the field 



right half space without the surfaces T, a . For Neu- </>(x) in the exponential of Eq. (11.10 1 is replaced by 
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d n (p(x) and G>(x,x',k) in Eq. (11.11 1 is replaced by 
(9„(9„'G>(x,x', k), where d n is the normal derivative 
pointing out of the objects. Here G>(x,x', k) is the free 
Green's function in the right half space which is given by 

G > (x,x', K )=G (x,x', K )tG (x,x^), (11.12) 

where x/j = (x||, —z) is the mirror image of x = (xm,z) 
and the — (+) sign applies to Dirichlet (Neumann) bound- 
ary conditions at the mirror plane at z = 0. The Green's 
function of free unbounded space is given by 



G (x,x', k) 



e ifc|x-x'| 

47r|x — x' 



(11.13) 



| k— >ik 



Using Eq. ( |H.12| the action S[g] defined by Eq. ( jll.llh 

can be expressed in terms of the original sources £> Q (x) 
and their mirror sources p^(x), 



S[g] = 




Ctx 1 Q* a (x) G (x, x' , K) Q/3 (x') 



J dxy , ^dx'^(x)G (x,x', K )gf(x') + c.c. | (II. 



14) 




FIG. 1: Geometry with original surfaces E a (objects D a ) in 
the right half space (z > 0) and mirror surfaces E„ (objects 
D^). x a , x^ are local coordinate vectors measured from the 
object's origins, Lp is the center-to- mirror distance of surface 
E/3 (object Dp). 



for Dirichlet boundary conditions at the surfaces E a and 
with Go(x,x',k) replaced by 9 n <9„<Go(x, x', k) for Neu- 
mann boundary conditions at the surfaces S Q . Here we 
have introduced the mirror sources 



0a ( x ) = Tea(xij), 



where the — (+) sign applies to a Dirichlet (Neumann) 
mirror. The mirror sources are located on the mirror sur- 
faces Ti^ that are obtained from the S Q by z — > — z for 
all surface positions, see Fig. [I] The first term of the ac- 
tion of Eq. ( |II.14 1 describes the interaction of the surface 
sources in the absence of the mirror plane. The second 
term couples each surface source to all mirror sources. 
Since the mirror problem is now described by an action 
in free space with sources and mirror sources, we can ap- 
ply the concepts of the previously developed approach 
for Casimir interactions between compact objects in un- 
bounded space (2] . Below, we provide an explicit deriva- 
tion for the case of Dirichlet boundary conditions at the 
surfaces S Q but we shall also indicate how the derivation 
has to be modified for Neumann boundary conditions. 
The action of Eq. (11.14 1 is composed of two quali- 



tatively different terms that we will now consider sepa- 
rately. Firstly, there are terms that couple sources on 
different surfaces (x ^ x') where we use the term "sur- 
face" in the following for the original and the mirror sur- 
faces. We shall call these terms off- diagonal. As diagonal 
terms we shall denote those which couple sources on the 
same surface (x — x' possible) . Both type of terms can 
be expressed in terms of the multipole moments of the 
sources. 



Off- diagonal terms — These terms couple sources on 
different objects, 



S a /3 — 



1 



dx. a {g* a (x a )</) /3 (x ct ) + c.c. } (11.16) 



(11.15) f or a ^ (3 and the original sources to the mirror sources, 



5'. 



a/3 



dx a {Q* a (x a )<pf(x a ) + C.C. } (11.17) 



for all a, [3. Here we have introduced local coordinates 
x a that are measured relative to an arbitrarily chosen 
origin that is located inside the surface a, see Fig. [T] We 
have also defined the fields 



18) 



x) = / dx'Go(x,x / ,K)e /3 (x / ) 
0f(x) = £^x'Go(x,x', K )^(x'), (II. 



which are the classical fields generated by the sources. 
For Neumann boundary conditions on the surfaces S Q , 



</>/3(x a ) and ^?(x„) in Eqs. (II.16I and (11.17) have to 



be replaced by d n c/)fj(x a ) and 3„rf(x a ), respectively. 



Also, Go(x,x',k) in Eq. (II.18I has to be replaced by 
<9„'G (x,x',k). 

Since we can assume that every position on S Q is out- 
side a sphere enclosing or E^, we can use the partial 
wave expansion of Go(x,x',k) for |x'| < |x|, 

G (x,x',k) = ji( iKr l h i 1 \ iKr ) Y hn(x)Y l * m (yi / ) . 



Irn 



(11.19) 
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When we consider Eq. (11.18) with coordinates x^ rel- 
ative to the origin inside surface Yip, the field that is 
generated by the source gp can be written as 

<j>p(xp) = -K^Qpj m hi(iKrp)Yi m (-ki3) , (11.20) 

lm 

where we have defined the multipole moments of the 
source gp as 



Qp,im= / dxpjt(iKrp)Y* m (xp)gp(-xp) . (11.21) 

The field of the mirror sources can be expressed in the 
same form, 



due to Eq. fllL15t and Y lm (k R ) = (-l)'+ m ^ m (x). For 
Neumann boundary conditions at the surfaces Y a the 
above expressions remain valid if the multipole moments 
Qp,im are replaced by Neumann multipoles which have 
the form of Eq. ( |IL21 1 but with ji{inrp)Y l * m {xp) replaced 
by d n \ji(iK,rp)Y£jxp)]. 



In order to express the action of Eq. (11.16 



in terms of 

multipole moments, we have to write the field generated 
by the surface source gp as a function of the local coordi- 
nate x Q that is regular at the origin inside Y a . This can 
be done using translation matrices U which relate outgo- 
ing {h^) and regular (ji) spherical Bessel functions, 



(xf) = -K^Q^h^iiKr^Y^if) , (11.22) 



lm 



where x^ denotes the local coordinates that are measured 
relative to the origin inside surface Y^. The multipole 
moments of the mirror source are given by 



Cp,lm 



T (-l) l + m Qp, 



(11.23) 



h\ \iKrp)Yi m (xp) = Ui> m 'im{'X- a p)ji> {inr a )Y Vm , (x a ) . 

I'm' 

(11.24) 

The matrix elements Ui' m 'i m of U depend on the vector 
X Q/ 3 = X^ — X a = x a — Tip from the origin inside Y a to 
the origin inside T,p, see Fig. [I] and are given by [2] 



Uv 



I'm'lrn 



(x 0/s ) = V^(-i) m i 1 - 1 ' V¥ 

x J^i 1 " ^21" + 1 



1)(2Z' + 1) 
I I' l"\ ( I V I" 

/ 1 772 —777,' 777' — 777 



h^){in\K a p\)Y v 



,(X a p), (11.25) 



where we have assumed that the Cartesian coordinate 
frames associated with the two origins have identical ori- 
entation, i.e., they are related by a translation. The 
summation over I" involves only a finite number of terms 
since the 3-j symbols vanish for I" > l + l' and I" < \l — l'\. 

Using the translation formula of Eq. (11.24 1, the 
field generated by the source on surface Yp, given by 



Eq. (11.20), can be written as function of the coordinate 
<pp(x a ) = -K^2Qp Mn ^2u^ llm ji>(iKr a )Yi, m ,(x a ) , 



lm 



where we have introduced W"^/ im = Ui 
ilarly we obtain for the field generated by the mirror 
sources of Eq. ( |II.22| , now expressed as function of x Q , 



(11.26) 
(X a( g). Sim- 



!>f(x Q ) = ±K^Qp Mn ^U^P lm jv{inr a )Y Vm ,{± a ) , 

(11.27) 



lm 



I'm' 



where we defined 



jjR,af3 1 1 \l+m 



I'm'lrn 



— (~ 1) m Ul'm'lm(x a 



4) 



(11.28) 



and used Eq. (11.23 1. Notice that the latter formula ap- 
plies also to the case a = /3 which describes the trans- 
lation between the surface Y a and its mirror image so 
that the argument of the translation matrix becomes 
x a — x^ = — 2L a z where L a is the normal distance be- 
tween the origin of surface Y a and the mirror plane. For 
this case of translations along the z-direction the trans- 
lation matrix simplifies to [2] 
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,(-2L Q z) = 5 m , m {-l) m i 1 - 1 ' ^{21 + 1)(2/' + 1) r'"(2i" + 1) 



I I' I" 




I I' I" 



-in 



Q 1 h%'{iK2L a ). (11.29) 



When we substitute the result for the fields of 



Eqs. (11.261, (11.271 into Eqs. (II.16I, (11.17 1, we obtain 



the action in terms of the original multipolc moments, 



Y = -« Y Y Qad'm' ^Vm'lm Qp 

a 7^0 Iml'm' 



I in 



(11.30) 



Yj^Sp — ± K Y< Y, Q*<x,Vm> tivm'lmQp.lm- 



a/3 



ct(3 Iml'm' 



(11.31) 



which applies to Dirichlet as well as Neumann boundary 
conditions at the surfaces E Q . To simplify notation, we 
have defined the modified translation matrix elements 



U alS 

I'm'lm 
U Vmlm 

I'm'lm 



(-1)%, 



Vmlr, 



(-if^L+c- 1 )^ 



Iml'm' 



(11.32) 
(11.33) 

(a ^(3) 
(11.34) 



Here we have used the definition of Eq. (11.21 1 and 
[ji(iKr)]* = (—\) l ji(inr). In addition, we applied 
the symmetry relations U l 



a/3' 
Irnl'r 



and 



Ut^l = {-l^'U^m which follow from Eqs. fIL25| 



(11.29 1, symmetry properties of 3-j symbols and spherical 

l ^2/(nz)K l+1/2 (z) where 



harmonics, and h 



n (jz) = —i -y £/ynz)I\ l + 1/2 

Ki + i/ 2 (z) is the (real valued) modified Bessel function 
of second kind. Notice that the actions of Eqs. (11.30 1, 



11.31 1 couple sources on different surfaces and hide the 



dependence on the particular boundary conditions and 
shape of the surfaces. 

Diagonal terms — These are the self-action terms 



Six 



1 



dx{p* (x)0 a (x) +c.c.} 



(11.35) 



in Eq. (II. 14 1, which we have expressed here in terms 



of the classical field defined in Eq. (II.18I. For Neumann 



boundary conditions on the surface S Q we have to replace 



again </> Q (x) by <9„0 Q (x) in Eq. (11.35 1. Here and in the 



following we only use the coordinate system associated 
with the origin inside S Q , and hence drop the label a 
on the coordinates. The classical field generated by the 
source g a (x) on surface E Q obeys the Hclmholtz equation 



(V 2 -K 2 )(/>«(x) = -g a (x). 



(11.36) 



For positions x that are located on the surface E a the 
total field </>(x) generated by all sources must obey the 



same Hclmholtz equation. The part of the total field that 
is generated by sources other than g a (x) can be regarded 
as incident field </)o, Q (x) at the surface E a , which obeys 
in region around S a , that is free of sources other than 
go,(x), the homogeneous Hclmholtz equation. Hence the 
total field can be written as 



(x) = 0o,o! (x) + (f) a (x) 



(11.37) 



= </> , Q (x)+ / dx'G (x,x',K)g Q (x') 

for all x located on S Q . For Neumann boundary con- 
ditions at S Q the Green's function in Eq. pi. 37 1 is re- 
placed again by d n >Go(x, x', k). We would like to eval- 



uate the action of Eq. (11.35 1 in terms of multipole mo- 



ments. Hence, we must consider field configurations with 
a fixed source on surface T, a that is characterized by its 
multipole moments. This implies that we have to find 
the incident field 0o,a(x) that induces a prescribed set of 
multipolcs Q a ,irn on E a . The multipole moments can be 
identified as the amplitudes of the scattered field which 
is given by Q (x) with x located outside of the surface 
E Q . Using the partial wave expansion of Eq. (11.19), we 
get 

0«(x) = -k Q a ,i'm'h'iP (inr)Y Vml (x) . (11.38) 



From scattering theory we know that the amplitudes of 
the scattered field are related to the amplitudes of the 
regular incident field by the transition matrix T which is 
related to the scattering matrix § by T = (§ — l)/2. If 
we expand the incident field as 

</>o,a(x) = Y ^o.cMm ji{inr)Yi m (±) , (11.39) 

tin 

the amplitudes of the scattered field are given by 

- KQa.I'm' = Y, T-Vm'lm </>0,a,im , (11.40) 



Im 



where the Ty m n m denote the matrix elements of the T- 
matrix of the surface E Q . Hence, the amplitudes of the 
incident field have to be given by 



t>0,a,lm - - K Y\T a }lml' m 'Qa 



l f m' 



(11.41) 



I'm* 



For Dirich let boundary conditions, the total field 0(x) of 
Eq. (11.37 1 has to vanish on S Q so that a ( x ) = — <?Va( x ) 



on the surface. Hence, using Eqs. (11.39 1 and (11.41 1, the 
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action of Eq. (11.35 1 can be expressed in terms of the 

(11.42) 



multipole moments of g a (x), 



Sat — K E Qa.lm, IX"] 



ImV 



Irnl'rn' 

with the modified T-matrix defined by 



j-ol _ / _i \ln-OL 
I'm'lm — V / I'm' 



Im ' 



(11.43) 



Here we have used the definition of the multipoles in 



Eq. (11.21 I to integrate over the surface and applied the 



relation T« m , lm = (-1)'+'' T^,. 

For Neumann boundary conditions on Yl a , the normal 
derivative of the total field of Eq. (11.37) has to vanish on 
E Q so that d n (f) a (x) = — <9 n <^ ,a( x ) on the surface. When 
we use the definition of the multipole moments for Neu- 



mann boundary conditions, we obtain again Eq. (11.42) 



but with the matrix T a for Neumann boundary condi 
tions. 

Now we can combine all results for off-diagonal and 
diagonal terms and express the total action of Eq. (11.14) 



in terms of the original multipole moments. Since we 
have 



a/3 



E*° 



E 

a,0 



(11.44) 



action of the multipole moments 



we obtain from Eqs. (|II.30|), (|II.31| and (]II.42j) the total 

(11.45) 



a/3 



Qf3 ■ 



where we have suppressed the sum over the indices I, to, 
I' to' and defined the matrix 



V af3 



k {[* T^a/s ~ U Q/3 (1 - M ± U^ Q/3 } (IL46) 



for Dirichlet (+) or Neumann (— ) boundary conditions at 
the mirror plane. The partition function of Eq. (11.11 1 is 



then obtained by integrating over all multipole moments, 



Z{k) = Z Yl / dQ a dQ* a exp 

a=l J 



Q/3 



na/3 



(11.47) 

The Gaussian integral over the multipoles is proportional 
to the inverse determinant of M Q/3 . Finally, we substitute 
into Eq. (II. 9 1 to obtain the Casimir energy, 



he 

2^ 



G?Kln - — — 
detfy 



(IE 



where the determinant is taken with respect to the par- 
tial wave indices I, m and the surface indices a, (3. The 
matrix is the result of moving the surfaces to infinite 
separation, where the translation matrices U, IJ R vanish 
so that 



JQ/3 



Q/3 ■ 



(11.49) 



In the special case of one compact surface in front of the 



mirror plane Eq. (11.48 1 simplifies to 



fir r°° 

Sx= — / d K lndet(l±T 1 U^ 11 ) 



(11.50) 



for Dirichlet (+) or Neumann (-) boundary conditions at 
the mirror plane. This expression applies to Dirichlet, 
Neumann and even more general boundary conditions 
at the compact surface which enter only via the matrix 
T 1 . Notice that this result depends on the original ma- 
trix elements (without tilde) since the phase factors of 
Eqs. (11.33 1 and (11.43) drop out when taking the matrix 



product of T 1 and U^' 11 . This general result shows that 
the Casimir interaction between a mirror and an object 
with arbitrary shape and boundary condition can be ob- 
tained from the transition matrix T 1 of the object and 
the translation matrix U^' 11 that describes the (classical) 
interaction between the induced source and its mirror im- 
age. 



B. Electromagnetic field 

The derivation of the Casimir energy for a scalar field 
can be extended to electromagnetic field fluctuations in 
the presence of dielectric objects pQ. The result will 
have the form of an effective action for electric and mag- 
netic multipoles of the current densities J Q inside the 
objects. We consider again N objects that are located 
in the right half space that is bounded by a perfectly 
conducting plane at z = 0. At his plane the tangen- 
tial electric field and the normal magnetic field vanish, 
E|| = 0, Bi = 0. The material objects are assumed to be 
dielectrics that are characterized by a frequency depen- 
dent dielectric function e a (u) and permeability function 

The partition function can be factorized again into a 
product of factors Z(k) at a fixed Wick rotated frequency 
k. Hence, in the following we consider all expressions at 
fixed k and suppress the label k. The Euclidean action for 
the electromagnetic field in the presence of macroscopic 
media without external sources can be expressed as 



1 



S = - dx [ED* + BH* + c.c. ] 



(11.51) 



in terms of the macroscopic fields D = e(in)E, H = 
/i _1 (zK)B. Here the energy density of the field is inte- 
grated since under a Wick rotation to imaginary time 
the Lagrangian in real time is generally transformed to 
the Hamiltonian in imaginary time. In this description 
the induced (bound) currents inside the material objects 
have been absorbed into the definition of the macroscopic 
fields. The partition function for this action is given by 
a functional integral over the vector potential A and the 
scalar potential <& (after introducing a Faddeev-Popov 
gauge fixing term) where the fields are expressed in terms 
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of the potentials as 



E = ikA-V$, B = VxA. 



(11.52) 



An alternative description in terms of the fields E, B 
only is obtained if the bound charges (p) and currents 
(J) density inside the objects are not substituted by the 
macroscopic fields but considered explicitly. Then the 
Wick rotated action can be written in terms of the po- 
tentials as 



S[3,p] 



dx 



* 2 (|A| 



|$| 2 ) + V|V^ 



(AJ* + <I>p* + c.c. ) 



IV*I 



(11.53) 



where we have chosen the Feynman gauge. The partition 
function is obtained by integrating over both potentials 
A, $ and sources 3, p. However, in the latter integration 



the currents and charges must be weighted according to 
the energy cost for inducing them on the objects. This 
energy cost must depend on shape and material of the 
objects. We will see below that this can be achieved by 
rewriting the self-energies of the separate objects in terms 
of the incident field that generates the polarizations and 
magnetizations which give rise to the induced current. 
Hence the action of Eq. ( II.53[ ) is independent of material 
and shape of the objects and these properties enter the 
partition function through proper weights on the currents 
that measure the susceptibility of the objects to current 
fluctuations. 



We proceed by integrating out the 
fluctuations of the potentials A, <f> 
Eq. pL53| 



unconstrained 
in the action of 
This integration yields the partition func- 
tion as a weighted functional integral over sources which 
we indicate at this stage by a subscript w on the integra- 
tion variable. When we denote the current density in the 
interior D a of object a by J Q , we get 



Z(k) = Z f[ f[D3 a ] w [D3* a } w ex P -^V / dx / dx' {3* a (x)G > (x 1 x' , k)3^) + c.c. } 

a =i J [ 2 a f> J Da Jd 

N r 

= [D3 a ] w [V3* a ] w e- s W, 



a=l 



(11.54) 



where we have used the continuity equation VJ = — inp to eliminate the charge density by introducing the tensor 
Green's function for the half space 



'10 

e>(x,x', K ) = Go(x,x / !K )I-G (x ) x / fl , K ) | 1 

,0 -1, 



1 



V ® V [G (x, x', k) - G (x, x' fl , «)] , (11.55) 



r 



where Gp( x, x', k) is the free scalar Green's function of 
Eq. (11.131, I is the identity matrix and x^ is defined 



below Eq. (11.12 I 



The action of Eq. 



11.54 1 can be expressed in terms of 



the original current densities 3 a and the mirror current 
densities. Their components parallel and perpendicular 
to the mirror plane are 

j£u(x) = -J Q ,||(x fl ), J« z (x) = J a , z {x R ) ■ (IL56) 

The action then reads 



dx / dx'J* (x)^ (x,x',k)J (3 (x') 




§[3] = 

+ J dxJdx.'3* a (x)g (x,x.',K)3f(x') + c.c.^ . 

(11.57) 

Now the sources are coupled by the free, infinite space 



Green's function 

£ (x, x', k) = G (x, x', k) + «) V'G (x, x', k) . 

(11.58) 

The mirror currents are located on the mirror objects 
that are obtained from the original objects by reflection 
at the mirror plane. This action has the same structure 
as in the case of scalar fields. It can be expressed in 
terms of multipolc moments of the current densities very 
similarly to the scalar case. Again, we consider diagonal 
and off-diagonal terms separately. 

Off- diagonal terms — We introduce again local coor- 
dinates x a that are measured relative to an origin inside 
D a . The terms that couple the original sources on dif- 
ferent objects can be written as 



S al 3 = w~ dx a 3* a (x a )E (x a ) +c.c. (11.59) 

J D„ 
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for a 7^ [3 and the terms involving mirror sources become which is an expansion in the outgoing solutions 



D af3 



1 



dx Q j;(x Q )Ef(x Q )+c.c. (11.60) 



2in j D 

for all a, /3. Here we have introduced the electric fields 
E^(x) = iK I dx'g (x,x',/t)J (3 (x') (11.61) 

J Dp 

Ef(x) = in j dx'0 o (x ) x',K)jf(x / ) (11.62) 
Jd* 

which arc generated by the current densities. For posi- 
tions x^ that are located outside a sphere that encloses 
the object Dp, the electric field can be expressed in terms 
of the electric and magnetic multipole moments of the 
current density Jp for I > 1, \m\ < I, 

Qiim = JdxpJp(x a )(-l) l Mj^(xp), (11.63) 
Qllm = |rfx /3 J /3 (x Q )(-l) i+1 N^(x /3 ), (11.64) 

where Mj m (x Q ), N; m (x Q ) are the regular, divergence- 
less solutions 

Mf m (x) = x [xj,(i/sr)Y, m (*)] (11.65) 

Nf m (x) = |-VxVx [xi ; (mr)y im (x)] (11.66) 

A IK 



of the vector Helmholtz equation and A = + 1). The 
electric fields of the original sources can then be written 
as 



E 



p(*p) = -iJ [Qilm M Uxp) + Q^ m Nf m ( X/3 ) 
lm 



(11.67) 



ML(x)-^Vx[x^ (1) (z K r)F im (x)] 



(11.68) 



1 1 



N? m (x) = ^--V x V x [x/i ; (1) ( iK r)r im (x)] (11.69) 

A IK 



of the vector Helmholtz equation. The electric fields of 
the mirror sources can be also expressed in terms of the 
multipole moments of the original sources, 



Ef(xf) 



lm 



[<f m ML(x£) + Q^>L(x£) 



(11.70) 

where x^ denotes the local coordinates of the mirror ob- 



ject D§. Using Eq. (11.56 1, the definition of the vector 



solutions of Eq. (|lL65| and Y lm (x R ) = (-l) l+m Y lm (x), 



the moments of the mirror currents are given by 



(-1) 



(11.71) 



Ql:L = -(-i) (+m < im - (n.72) 

The expansion of the electric fields into outgoing vec- 
tor waves with respect to the origin of the object that 
generates the field does not allow us to perform the in- 
tegrations in Eqs. (11.59), ( II.60[ ). We would like to ex- 
pand the electric field generated by object Dp in terms 
of vector waves that are regular at the origin of object 
D a . When the coordinate systems associated with the 
two objects have identical orientation, this can be done 
by relating the outgoing and regular vector waves by a 
translation matrix U. Generalizing the result for scalar 
fields, the translation matrix couples both types of vector 
solutions, 



M 



L(x/3) = ^{ J B rm ^ m (X Q/3 )Mf, m ,(x Q ) + Cp m - ;m (X C[/3 )N^ m ,(x Q )} 

I'm' 

NL(x/3) = {Sl'm'Zm(X Q/ 3)N£ m ,(x Q ) + Q' m ' lm (X a p)Mj, m , (* a )} 



(11.73) 
(11.74) 



I'm' 
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with X Q) g =x„ — X0. The matrix elements are given by |35j . 



;l" 



'(2Z+1)(2Z' + 1)(2Z" + 1) 
/(/ + !)/'(/' + !) 



i r i"\ f i v i" 

/ \m —m! ml — m 



hy) (iK\X. a p\)Y v , m - m i (X Q/3 ) 



(11.75) 



^2 (^toi^ i ' m ' i » ra + 1 (-^- Q / 3 ) + ^; m J 4i'm'im-l(X Q ^)) 



im A 



V m'lm 



( X a/3) 



(11.76) 



with 



A, Wi m(X Q/9 ) = (-l) m i'-'' ^ ^(2/ + 1)(2Z' + 1)(2Z" + 1) 

l" 

X (o o) (m -m' r/ m) ^(-IX^D^^C^) (0.77) 



and = ^(i^mjflim+l). It is useful to combine 
the translation matrix elements to the matrix 



7J a/3 = 



that acts on the vector Qf m = (Q° lm , . 
moments. 



(X Q/3 ) (11.78) 



Qe im) of multipole 



Using the translation formula of Eqs. (11.73), (11.74 1 



we can express the field E^(x Q ) in Eq. (11.59) in terms 
of regular vector solutions and perform the integration. 
This integration yields the multipole moments defined by 
Eqs. (IIL64T), dlL63T). The action of Eq. pL59l can then 



be written as 



E 



=E E 

Iml'm' 



Vm' U V 



I'm'^l'm'lm^lm 



Q 



(11.79) 



where we have defined the modified translation matrix 



I'm'lm 



i 7 y/3« 



(11.80) 



where f denotes the conjugate transpose of the ma- 
trix of Eq. (11.78). For translations along the z-axis, 
~K. a f3 ~ z, the latter expression is diagonal in 77i 
and simplifies to Uf mlm = {-^) l 'uf mlm due to the 
symmetries Bf ml , m (Xf3 a ) = (-l)' +r B rmim (X Q/3 ) and 

C *TOi'm( X /3<*) = -(-1)' +/ CVmim(X a/3 ). 



The action of Eq. (11.60) that couples the original 
and mirror sources can be cast into a form similar to 



Eq. (11.79 1 with a translation matrix V R that is defined 

by 



j,R,al3 
"I'm'lm 



= (-1) 



l-\-m 



-^I'm'lm 
~C[> 'm'lm 



-C v 



m'lm 



(11.81) 

where the phase factors relating the original and mirror 



multipole moments in Eqs. (11.71 1, (11.72) have been ab- 



sorbed in the translation matrix. 
Eq. ( |O.70| becomes 



With this definition 



E^ 

a/3 



^ QMtlLQL (n.82) 



IrnV m' 



with the modified translation matrix 



fjR.a/3 
I'm'lm 



1 r 

2 



(-i) ; «L+(-im 



Iml'm' 



(11.83) 



For a = P, the action of Eq. (11.82 1 describes the interac 



tion between a source and its mirror image. In this case 
the translation vector is x a — = — 2L a z where L a is 
the normal distance between the origin of object D a and 
the mirror plane. Hence, the translation is along the z- 
axis and the translation matrix elements of Eqs. (11.751, 
(11.76 ) simplify to 
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-I" 



B Vm , lm (-2L a i) = 5 m , m {-l) m i 1 - 1 ' ]T V + 1} + + + 1)] 



C/'m'im(~2Z/ Q z) 



I I 1 l"\ (I VI 



/ \ m —m 



oJU -m oJW" 22 ^ 

I 



'(2Z + 1)(2Z' + 1) 



(2/" + 1) 



(11.84) 



(11.85) 



These matrix elements obey the symmetry relations 
B i> m i m (- 2L az) = Bi mVm {-2L a i), Cp mlm (-2Lg Z) = 
Ci m i' m (—2L a z) so that the matrix of Eq. ( |II.83 1 sim- 
plifies for a = (5 to W*X = (-1)%^°° 



I'mlm ' 



Diagonal terms — So far we have described the in- 
teraction between arbitrary multipolcs inside the mate- 
rial objects. But the functional integral over currents 
of Eq. (11.54 1 contains weights that measure the energy 



cost for inducing currents on an object with particular 
shape and material composition. In the following we will 
show that these weights can be implemented straightfor- 
wardly when we express the diagonal terms of the action 
of Eq. (|IL57|, 



S a — 



c.c. } 



I [ dx [ dx'{J* t (x)£ (x,x', K )J ce (x') 

(11.86) 

in terms of the incident field which generates a prescribed 
current that corresponds to the induced polarization and 
magnetization inside the object. Then the relation be- 
tween the induced current and the incident field ensures 
that the currents arc weighted properly. To see this we 
employ the macroscopic formulation of the electromag- 
netic action, see Eq. (11.51 1. The energy of Eq. (11.86) 



is associated with the process of building up the induced 
current J Q inside the object. This energy must there- 
fore be equal to the change in the total macroscopic field 
energy that results when the object is placed into an ex- 
ternal field. Thus we obtain 



§ 1 
ba ~ 2 



rfx [ED* + BH* - (E Dq + B Hq) + c.c. ] 



(11.87) 

where integration extends over all space. Field vectors 
with subscript represent the incident field that is gen- 
erated by some fixed external sources in otherwise empty 
space and field vectors without label stand for the total 
field from the external and induced sources after adding 



the object. S a can be also written as 
1 



S a = -ydx[ED3-D*E + BH5-H*B + c.c.] 

+ \ |dx[(E + E )(D*-D;) 

+(B + B )(H* - HS) + c.c. ] . (11.88) 

The second integral of this expression vanishes. This can 
be seen by setting E + E = inA - V<I>, B + B = V x A 
and using that VD* = p* xt = VD*, V x H* + inT>* = 
Jcxt = V x Hq + ikDq since the external sources p ex t, 
J cx t are fixed when the object is added. Notice that the 
bound (induced) sources do not appear explicitly since 
they are included in the macroscopic fields. Since D = 
E , H = B and D = eE, H = ^i~ 1 B with e = e a , 
/.t = fi a inside the objects and e = 1, fj, = 1 outside the 
objects, we get 



/ dx[(l-e Q )EE; + (l-/i- 1 )BBS+c.c. 



s 1 

2 

(11.89) 

where integration runs only over the interior of the ob- 
ject. The material dependent functions e a and [i a can 
be expressed in terms of the polarization P a and mag- 
netization M a of the object. The relations between 
macroscopic fields yield P a = D E = (e a — 1)E and 
M a = B - H = (1 - l//x a )B. When we substitute 

to obtain 



i/nV x EJ in Eq. ( II.89[ ) we can integrate by parts 



<7x 



-J a (x)Eo(x) + c.c. 



(11.90) 



where we have combined the polarization and magneti- 
zation to yield the induced current density 



V x M Q . 



(11.91) 



Notice that the incident field E in Eq. (11.90 1 depends 



on the current density 3 a since Eo has to induce the pre- 
scribed current density. Hence, the problem of expressing 
the diagonal part of the action in terms of multipole mo- 
ments has been reduced to computing the incident field 
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that has to act on the object to generate a given current 
density. In scattering theory one usually encounters the 
opposite problem. For an incident field one would like 
to compute the scattered field which can be expanded 
in outgoing partial vector waves, see Eq. (11.67 1. Here 
the situation is slightly different. We seek to determine 
the incident field that generates a given set of multipole 
moments inside the object. In other words, for a given 
scattered field, which according to Eq. (11.67) is given by 



the multipole moments, we would like to obtain the cor- 
responding incident field. We expand the incident field 
as 



E (x) 



^ 2 E 

I in 



(11.92) 

The relation between the multipole moments and the 
amplitudes of the incident field is determined by the T- 
matrix T° = (§ Q — 1)/2, where S" is the scattering matrix 
of the object. When we solve this relation for the incident 
field amplitudes <j> lm = (4> M ,im, ?Wm) 5 we get 



4>lm — X]P~ a ]imZ'm'Q"m' ' 



(11.93) 



I'm' 



where Tl m ii m i is a 2 x 2 matrix acting on magnetic 
and electric multipoles. This relation together with 
Eq. (11.92 1 yields the incident field that generates the 



multipoles Q^ n . With this result we can perform the in- 
tegration in Eq. (11.90 1 which yields S a in terms of the 
multipoles, 



F 

Iml' 



i 



[T 



on-i 

Iml'm' 



[fa] 



I'm'lm 



with the modified T-matrix 



-T, 



1 lmV 



(-1)' 



Iml'm' 



Iml'm' 



T, 



Iml'm' 



rpOt,EE 

Iml'm' 



Ql'm' 
(11.94) 

(11.95) 



where the Tf*^ , are the elements of the matrix T^,^ 
with X, Y — M, E labeling magnetic and electric ele- 
ments, respectively. Due to the symmetry 



//TICK, MM rpOt, ME 

Iml'm' Iml'm' 



1 Iml'm' 1 I 



("I) 



l+V 



rpOt, MM 

I'm'lm 



mQ , EM 

I'm'lm 



Iml' 



-T° 



,'lm 



rr,a,EE 

I'm'lm 



of the T-matrix Eq. ( [II.94 1 can be simplified to 



Iml'm' 



(11.96) 



(11.97) 



The total action S[Q] for the multipole moments is 
given by the sum of the actions of Eqs. (11.79 1, (11.82) 
and (11.97 1. In compact matrix notation it can be written 
as 



5[Q] = FQ« r 

a/3 



ia/3 



Q 



p. 



(11.98) 



with the matrix 



where T° 



(11.99) 

\J a P and \J R ' a P stand for the matrices with 



matrix elements given by Eqs. (11.95 1, (11.80 ) and 
respectively. 



11.83 



11.47 



In analogy to the scalar case, Eq 
the partition function Z(k) is obtained by integrating 
over all multipoles. Notice that by construction of the 
diagonal parts of the action of Eq. (11.98 1 the proper 



weights are assigned to the multipole configurations. 



The Gaussian integral over multipoles and Eq. (II. 9 1 
lead to the final result 



£ = 



he 

2k 



duln 



detl 
dctIV 



(11.100) 



for the electromagnetic Casimir energy where the deter- 
minant is taken with respect to the partial wave indices 
I, m, polarization indices M, E and the object indices 
a, 0. The matrix is given by Eq. dlL99| with the 



translation matrices tS a P and U^'"^ set to zero. 

For one object in front of a perfectly reflecting mirror 
plane, the Casimir energy can be written as 



Si 



he 

2n 



d/t bidet (l + fMD*' 11 



(11.101) 



The matrix U^' 11 that describes the interaction between 
fluctuating currents on the object and its mirror image 
is universal, and the shape and material composition of 
the object enters through its T-matrix T 1 . 



III. PLATE - SPHERE GEOMETRY 

In this Section we consider a geometry that is most 
relevant to a large number of experimental studies of 
Casimir interactions carried out in the last decade. This 
geometry, consisting of a plane mirror and a sphere, is 
experimentally favorable since it avoids the problem of 
parallelism for plane surfaces facing each other. Despite 
its experimental importance, this geometry lacks a theo- 
retical description of the electromagnetic Casimir inter- 
action. For a scalar field with Dirichlet boundary condi- 
tions at the sphere and the plane, the interaction over 
a wide range of distances has been obtained recently 
[2U1 IBUj . including an analytic expression for the low- 
est order correction to the proximity force approximation 
[38]. 

We consider a sphere of radius R in front of a plane 
mirror with distance L between the center of the sphere 
and the mirror, see Fig. [2] For the scalar field we study 
Dirichlet or Neumann boundary conditions at the mirror 
and the sphere. For the electromagnetic field, the sphere 
consists of a material with arbitrary dielectric function 
t(k) and permeability /i(fc) while the mirror is assumed 
to be perfectly reflecting (implying vanishing of the par- 
allel components of the electric field and the normal com- 
ponent of the magnetic field) . 
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FIG. 2: Length scales of the sphere-mirror geometry. 



A. Large distance expansion 

We start by studying the Casimir interaction at large 
separations. For separations L that are large compared 
to the size R of the object, the Casimir energy can be 
expressed as an asymptotic series i n R/L. Using that 
lndet = tr In in Eqs. ( |II.50| | and ( fll.lOl) , we get by 
expanding the logarithm 



he 

2^ 



cfctr ln(l - N) 



(HI.l) 



he 

"2?r 



J dK Y, p tr ( w )> ( IIL2 ) 



where we have defined the operator N = ipT 1 !^' 11 with 
— (+) for a scalar field in the presence of a Dirichlet (Neu- 
mann) mirror and N = — T 1 !^' 11 for the electromagnetic 
field and a perfectly reflecting mirror. The operator N de- 
scribes a wave that travels from the mirror to the object 
and back, involving one scattering at the object. 



1. Scalar field 

All information about the shape of the surface and the 
boundary conditions is provided by the T-matrix. For a 



spherically symmetric surface the matrix is diagonal and 
completely specified by scattering phase shifts Si(k) that 
do not depend on to, 



Tlml'm'{k) = Su'S, 



;l< in' 2 ^ ' 



2k5, (k) 



- 1 



(III.3) 



where k is the real frequency. The scattering phase shifts 
for a sphere of radius R with Dirichlet and Neumann 
boundary conditions are 



cot Sf(k) 
cot 8f(k) 



ni(H) 

m 
m 



(III.4) 
(III.5) 



where £ = kR and ji (to) are spherical Bessel func- 
tions of first (second) kind. To make use of the gen- 
eral result of Eq. ( II.50[ ), we have to evaluate the matrix 
elements for imaginary frequencies k = in. This can 
be done using ji{iz) = i l y/ir /(2z)/, +1 / 2 (z), /i ; (1) (iz) = 

jl(iz) + ini(iz) = -i~ 1 v / 2/(irz)Ki +1 / 2 (z), where I/+i/ 2 
and K l+ i/ 2 are modified Bessel functions of first and sec- 
ond kind, respectively. We obtain 



Imlr. 



q-N 

Imlr. 



, 1 y7T h+l/l{z) 

1 ' 2 K l+1/2 (z) 

f 7T h+l/2{z)-2zl[ +1/2 {z) 
l + l/2\ Z > 



(III.6) 



2K l+1/2 (z)-2zK' H .a( IIL7 ) 



for Dirichlet and Neumann boundary conditions where 
z = kR. 

The distance dependence of the Casimir energy of 
Eq. (11.50 1 is given by the translation matrix that de- 



scribes the interaction between the original and mirror 
multipole moments. Since the origins (centers) of the 
original and mirror sphere are related by a translation 
along the z-axis, the matrix elements of the translation 
matrix U^ 11 are given by Eqs. pl.28|) and (|II.29|. After 



a rotation of frequency to the imaginary axis we obtain 



I'mlm 



-(-l)^- r +V(2^+l)(2/' + l) (-!)'" (2f" + 1) 



I I' I" 




/ /' I" 

771 —TO 



K v , +1/2 (2kL) 
V ttkL 



(III. 



r 



With the elements of the transition and translation ma- 
trices we can compute the Casimir interaction from 
Eq. (11.50 1. The scaling of the T-matrix for small k 



shows that partial waves of order / start to contribute 
to the energy at order L _2 ^ +1 ' if the T-matrix is diag- 
onal in /. Hence, we can truncate the infinite matrix 



in Eq. (III.ll at finite multipole order to obtain the 



large distance expansion of the energy at a given order in 
R/L. Also, we can truncate the series over p in Eq. ( III. 1 1 
since the p th power of N becomes important only at order 
The first p powers of N describe p scattering off 
the sphere. 

We have used Mathematica to perform the matrix op- 
erations and to expand the integrand of Eq. ( III. 1 ) in 
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R/L. From this we find that the Casimir energy can be 
written as 



J'=2 



(III.9) 



where bj is the coefficient of the term ~ . In the fol- 
lowing we give the results for a mirror plane with Dirich- 
let or Neumann boundary conditions separately. 

Dirichlet mirror — In this case we have N = — T 1 U' R ' 11 . 
If the field obeys Dirichlet boundary conditions at the 
sphere, we obtain the for leading coefficients 



b 2 = 
b 5 = 
b* = 



1 



3083041 



421 
1152 ' 



535 

9216' 6 ~ ~ 4147200' 
557222415727 



292626432000 



2741117 
11059200' 

(111.10) 



indicating an attractive force. This result agrees with 
recent findings of Wirzba |30j . If we impose Neumann 
boundary conditions at the sphere, we find 



b 2 = 0, 6 3 = 0, 6 4 = 
b s 



17 

96' 



= 0, 6 6 = 
243487 



379 
960' 



b 7 = 



479 
24576 ' 



403200 



(III.ll) 



This result corresponds to a repulsive force which is ex- 
pected for unlike boundary conditions. To leading order 
the energy scales as R 3 /L 4 since the first two coefficients 
vanish. This behavior can be understood from the ab- 
sence of low- frequency s-wave scattering for a sphere with 
Neumann boundary conditions. 

Neumann mirror - In this case we have N = 
+T 1 V R ' 11 . For a Dirichlet sphere we obtain the coef- 
ficients 



bi 
b 5 
b S 



1 



64 = 



403 



64' 1152' 

73 , 2981791 

, Or = , 07 

9216 4147200 

550781720977 



292626432000 ' 



3636227 
11059200' 

(111.12) 



and hence the force is again repulsive. The modulus of 



the coefficients are smaller than those in Eq. (III. 10 1 for a 
Dirichlet mirror with the exception of the leading one, 62 



and the one with inverted sign, 67. Notice that the sum 
of the Casimir energies for the Dirichlet and Neumann 
mirror, both opposite to a Dirichlet sphere, is identical 
to the Casimir energy between two Dirichlet spheres at 
a center-to-center distance 2L [29] . The reason for this 
is that any field configuration can be decomposed into 
symmetric and antisymmetric modes with respect to the 
mirror plane which obey Neumann and Dirichlet bound- 
ary conditions on that plane, respectively. Since the en- 
ergy between two spheres scales as R 2 /L 3 for asymptot- 
ically large separations, the modulus of the coefficients 
62 in Eqs. (III. 10 1 and (III. 12 1 must be identical so that 
the contribution ~ R/L z cancels. It is easily seen that 
the higher order coefficients of Eqs. (III. 10 1 and (III. 12 1 



combine to the correct coefficients of the large distance 
expansion for two Dirichlet spheres [2J. 
For a Neumann sphere we get 



b 2 
h 
h 



0, 64 

_379 
~960' 



17 

96' 



267 
8192' 



243487 
403200 



(III. 13) 



For the same reason as before, the first two coefficients 
vanish. Notice that all given coefficients, with the ex- 
ception of 67, are equal to minus the corresponding co- 
efficients for a Neumann sphere and a Dirichlet mir- 
ror in Eq. (III.ll I. This result is consistent with the 



observation that the Casimir energy for two Neumann 
spheres at a center-to-center distance 2L is given by 
the sum of the energies of a Neumann sphere at dis- 
tance L from a Dirichlet mirror and a Neumann mir- 
ror, respectively. Since two Neumann spheres interact at 
large distance with an energy ~ R 6 / L 7 with the next- 
to-leading term ~ R s /L 9 , the sum of the corresponding 
coefficients, excepting 67, must vanish. The coefficients 
67 in Eqs. (III.ll I, 13[ ) combine to the correct value 
for two Neumann spheres that was obtained in Ref. [2] . 



2. Electromagnetic field 

We consider a dielectric sphere at a center-to-surface 
distance L from a perfectly conducting mirror. Due to 
spherical symmetry, the electric and magnetic multipolcs 
for all Z, to are decoupled so that the T-matrix is diagonal 



11, 



rpMM 
Imlm 



(-i)'~ 









f 2t7z/; +; 


(nz) 


- nl l+ i(nz) 




t2zl' l+i (z) 






*(*) 


I l+ i(nz) ^ 




(nz) 


-nl l+ i(nz) 


K l+h (z) 


±2zK> +i (z) 



(III. 14) 
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where the sphere radius is R 



Tfmim is obtained from Eq. 



kR, n = \J e(iK)/i(iK), 
( III.14[ ) 

by interchanging e and [i. The limit of a perfectly con- 
ducting sphere is obtained by taking e(ift) — + oo at an 
arbitrarily fixed [i(in) which can also vanish. Then the 
matrix elements become independent of fi, 



Imlr. 



1 ] 2^W 



(III. 15) 



r, 



I mini 



I l+1 {z)+2zl' k {z) 

= (-1)'- 2 — ^ (III 16) 

1 J 2^ +i »+2z^' +i (z) 1 DJ 



For all partial waves, the leading low frequency contri- 
bution is determined by the static electric multipole po- 
larizability, of = [(e - l)/(e + (I + l)/l)]R 2l+1 , and the 



corresponding magnetic polarizability, af = [(/i — l)/(/i+ 
(l + 1) /l)]R 2l+1 . Including the next to leading terms, the 
T-matrix has the structure 



Imlr. 



.21 



(-iy-\l + l)af 
l{2l + l)V.(2l-l)\\' 



'7i4 K 



^d T™ lm 



is obtained by af 



a l > llr, 



it- The 

first terms are 7f 3 = -[4 + fi(efi + fj, - 6)]/[5(/i + 2) 2 ]i? 5 , 
7 f 4 = (4/9)[(/i-l)/(A*+2)] 2 i? 6 , and 7^3, jf 4 are obtained 
again by interchanging /i and e. Higher order terms can 
be easily obtained by expanding Eq. (III. 14 1 for small n. 
The 



translation matrix elements 



l' 7 ,R,U 



ILai 



tin 



are obtained from Eqs. (11.81), (11.84 1 



and (11.85). Using h t 
get 



(1: 



f.RM 
u l'rnlm 



(iz) = -i 1 v / 2/(ttz)K 1+1/2 (z) we 



Bl'm'lmi — 'ZLz) 



(-1) 



C Vm , lm {-2Li) = 5 m , m 2KLm{~l) m i 1 - 1 ' Y,i- l f ' 



I V V 



I I' l"\ K t 



/ \m — 771 



H" + l/2_ 

V ttkL 



, {2l + l){2V + 1) 



Z J' f"\ / J V l"\ K v , +1/2 (2kL) 
O) [m -m O) ' 

r /(2Z + l)(2/' + l) 



J(f+l)J'(f' + l) 
(2kL) 



(21" + 1) 



(2/" + 1) 



(III. 17) 



(III. 18) 



Now we can employ the series representation of the 
Casimir energy in Eq. fall.ll with N = -f 1 U R,U . Since 



the T-matrix is diagonal in /, partial waves of order I 
start to contribute to the energy at order L~ 2 ^ l+1 \ Also, 
the p th power of N becomes important only at order 
L~( 3p+1 >. Notice the stronger increase of the exponent 
with p compared to the scalar case where the exponent 



is — (p + 1). This can be understood from the absence of 
s-waves for the electromagnetic field so that each reflec- 
tion contributes a factor 1/L 3 due to p- waves. Matrix 
operations and the expansion in R/L are performed with 
Mathematica. From this we obtain for a dielectric sphere 
in front of perfectly conducting mirror plane the Casimir 
energy 



he { 3 



;K-<) Z 4 + 



1 



15 
32 
7 



"2 



1 



+ 27^-27-)-,+ 



1 



L 6 1024 



[23«) 5 



Uafofi 



+ 2160(714 " 7i M 4 )] Jv + Y^ t 572 ( a 3 " «3 ) + 675 (9( 7l E 5 - 7 f 5 ) 



55( 72 E 3 



7^3))] 



L 8 



E^2 



23K) 



(III. 19) 



The electric contribution to the leading term, ~ L~ 4 , 
was obtained by Casimir and Polder for the interaction 
of an atom with static polarizability a\ and a metal- 
lic surface [TS]. Later Boyer has generalized the lead- 
ing order result to include magnetic effects described by 
af |36j . The higher order terms are new. They show 



how higher order polarizabilities and frequency correc- 
tions to the static parameters influence the interaction. 
There is no ~ 1/L 5 term. Notice also that the first three 
terms of the contribution at order L~ 7 have precisely the 
structure of the Casimir-Polder interaction between two 
atoms with static dipole polarizabilities af and af but 



1G 



it is reduced by a factor of 1/2 8 . This factor and the 
distance dependence ~ L~ r of this term suggests that it 
arises from the interaction of the dipole fluctuations in- 
side the sphere with those inside its image at a distance 
2L. The additional coefficient of 1/2 in the reduction fac- 
tor (l/2)(l/2 7 ) can be traced back to the fact that the 
forces involved in bringing the dipole in from infinity act 
only on the dipole and not on its image |37j . 

For a perfectly conducting sphere the coefficients of the 
expansion of the Casimir energy in R/L become univer- 
sal numbers. These coefficients can be cither obtained 



from Eq. (III. 19 1 by using the appropriate values of the 



parameters for a prefect metal or it can be computed di 
rectly f rom the T-matrix elements given in Eqs. (III. 15 1 
(III. 16 1. Following the latter route, we get the series 



OC 



he l ^ 

L ^— ' 

J =4 



7T 



(111.20) 



where the coefficients up to order l/L 11 

9 h n h 25 
16' 65 = 66 = "32' 
12551 . 1282293 



64 = - 



b 7 = - 



9600 ' 69 
32027856257 



no 



163840 ' 

39492614653 

on - 



3023 
4096 



(111.21) 



722534400 412876800 
This and the corresponding results for a scalar field ap- 
pear to be asymptotic series. Therefore, the series cannot 
be summed to obtain the interaction at small separations. 
In the next Section we use a numerical implementation of 
Eq. (11.101 ) to compute the interaction at all separations. 



B. Non-perturbative result at all separations 

To obtain the Casimir interaction over a broad range of 
distances, Eq. (11.101 1 has to be evaluated numerically. 



As we employ a partial wave expansion, the computa- 
tional work increases with decreasing separation between 
the objects. However, we shall see below that even at 
small separations our method yields sufficient precision to 
obtain the leading corrections to the proximity force ap- 
proximation (PFA). The numerical approach is based on 
the technique presented in Refs. [HE]. Using the analytic 
expressions for the matrix elements of the translation and 
transition matrices, we compute the determinant and the 
integral over frequency k in Eq. (11.101 ) numerically The 



matrices are truncated at a finite multipole order I which 
yields a series of estimates £™ for the Casimir energy. 
The exact Casimir energy is then obtained by extrapo- 
lating the series to I — > 00. For the geometry considered 
here, we observe an exponentially fast convergence that 
allows us to obtain accurate results even at small sepa- 
rations from a moderate multipole order. 

We compare our results for the interaction energy 
to the estimate that follows from PFA. This is impor- 
tant since PFA is often used over a range of separa- 
tions although its accuracy is unknown even at small 




FIG. 3: Sphere-plane Casimir interaction of a scalar field with 
Dirichlet boundary conditions on both surface s. The energy 
is scaled by the PFA estimate of Eq. (|III.22b. The dashed 



curve represents the large distance expansion of Eq. (III. 9 1 
with the coefficients of Eq. ( |III.10[ ). Inset: Corrections to the 
PFA at small separations as function of the surface-to-surface 
distance d — L — R. The dashed curve corresponds to the 
lowest order correction to the PFA obtained by Bordag 38 . 



distances for most geometries, including the electromag- 
netic Casimir interaction between a plane and a sphere. 
The PFA estimate for the latter geometry is given by 



£pfa = 7T* 



±/EM 



hcR 



(m.22) 



where d — L — R is the surface-to-surface distance. The 
amplitudes follow from the result for two parallel plates 
and are given by 



*0 + 



1440 
7tt 2 
11520 

720 



(111.23) 
(111.24) 
(111.25) 



for a scalar field with like (— ) and unlike (+) boundary 
conditions and the electromagnetic field, respectively. In 
the following we present our results for the plane-sphere 
interaction over a wide range of separations for scalar 
fields with Dirichlet and Neumann boundary conditions 
and for the electromagnetic field. 



1. Scalar field 

We consider the four cases that correspond to Dirich- 
let or Neumann boundary conditions at the sphere and 
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FIG. 4: Analog of Fig. [3] for Neumann boundary conditions 
at the sphere and the plane. The dashed curve represents the 




large distance expansion of Eq. (III. 9 1 with the coefficients of 
Eq. dliTTl). 



FIG. 6: Analog of Fig.[3]for Dirichlet boundary conditions at 
the sphere and Neumann boundary conditions at the plane. 
The dashed curve represents the large distance expansion of 
Eq. 1III.9I) with the coefficients of Eq. (|III.12|). 
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FIG. 5: Analog of Fig. [3] for Neumann boundary conditions 
at the sphere and Dirichlet boundary conditions at the plane. 
The dashe d curve represents the large distan ce expansion of 
Eq. 1111.91 with the coefficients of Eq. dlll.lll. 



the mirror plane. The results for the Casimir energy are 
shown in Figs. 3][6 At large separations the results of the 
numerical evaluation of the determinant agree nicely with 
the large distance expansions of the previous section. For 
a sphere with Dirichlet boundary conditions, the energy 



scales both at large and small separations as R/d 2 so 
that the curves for £/£pfa in Figs. [3] and [6] tend also to 
a constant for R/L — > 0. The PFA underestimates the 
actual interaction energy at all separations. For a Neu- 
mann sphere the interaction scales at large separations 
as R 3 /L 4 so that the curves for £/£pfa in Figs. [3] and [S] 
tend to zero for R/L — > 0. Here the PFA overestimates 
the actual interaction energy at all separations. 

The non-perturbative approach allows us to study also 
the case of small separations. In that limit our results 
can be fitted to a power law of the form 



£ — £pFA 



R 



(111.26) 



The coefficients 6j obtained from a fit of the function 
of Eq. (III. 26 1 to the data points for the four smallest 
studied separations are summarized in Tab. [T] The fitted 
curves are shown as insets in Figs. |3][6| The result for 
6>i for the case of Dirichlet boundary conditions at the 
sphere and the mirror agrees with the analytical result 
B\ = 1/3 presented in Ref. [38]. Again for Dirichlet con- 
ditions, for the second order coefficient a considerably 
larger numerical estimate of 62 — 1.92 has been obtained 
from world- line Monte Carlo sampling [31] [32] . For the 
other combinations of boundary conditions our findings 
represent the first results for the corrections to the PFA. 
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boundary condition 


0i 


02 


plane / sphere 






D / D 


0.3346 ±0.0017 


-0.231 ±0.012 


N / N 


-2.4298 ±0.0070 


4.394 ±0.049 


D / N 


-1.7493 ±0.0061 


2.354 ±0.042 


N / D 


0.1124 ±0.0043 


1.240 ±0.044 


EM / EM 


-1.42 ±0.02 


2.39 ±0.14 



TABLE h Coefficients describing corrections to PFA, see 
Eq. ftll.26) . D=Dirichlet, N=Neumann (scalar field) and 
EM=perfectly conducting boundary conditions (electromag- 
netic field). 



2. Electromagnetic field 

We focus on a perfectly conducting sphere and mirror 
plane. The Casimir energy resulting from a numerical 
computation of the determinant of Eq. (11.101 1 is shown 
for a wide range of separations in Fig. |7J At large sep- 
arations the interaction is described by the asymptotic 
result of Eq. (III. 20 1. At small distances the Casimir en- 



ergy approaches to PFA estimate and corrections to PFA 
can be described again by Eq. (III. 26 1. A corresponding 



fit to the data points for the four smallest separations is 
shown as inset in Fig. [7] The corresponding amplitudes 
of the correction terms are listed in Tab. |U 




This result is important for a number of recent mea- 
surements of the Casimir force between almost flat sur- 
faces. These experiments have been performed for a 
plane mirror and a sphere with a radius that is much 
larger than the distance between the surfaces in order 
to avoid difficulties from parallelism control. Our results 
at short distances indicate that for a sphere of radius 
R = 100/im the corrections to PFA are below 1% only 
for surface-to-surface distances d that are smaller than 
sa 700nm. For a sphere of that size, PFA fails already by 
10% at a distance of » 8/xm. From this we conclude that 
the currently achieved experimental accuracy for mea- 
surements of Casimir forces in the sphere-plate geometry 
with a sphere of R = 100/im at distances below d = 1/im 
is within the range or slightly less than the corrections 
to PFA. 



However, deviations from PFA become severe when 
smaller objects interact with surfaces. Hence it is im- 
portant to understand the crossover of the Casimir in- 
teraction between macroscopic objects and the eventual 
Casimir-Polder interaction between single atoms and a 
surface. A description of this crossover is provided by 
our results in Fig. [7] Corresponding results can be also 
obtained for a dielectric sphere or less symmetric objects 
by the methods presented here. Finally, we note that 
it should be also possible to compute correction ampli- 
tudes like the 8j in Eq. (III. 26 1 analytically by applying 
methods similar to those used in Ref. |26j . 



However, the 

validity range of the corresponding series would be lim- 
ited to small separations and for an overall description 
of the interaction one should resort to the device of a 
numerical evaluation along the lines presented here. 



0.2 0.3 0.4 0.5 
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